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Abstract 

■ Exact linear retarding spherically symmetric solutions of Einstein equa- 

' tions linearized around Friedmann background for the ultrarelativistic 

^ , equation of state are obtained and investigated. Uniqueness of the so- 

I ^ ' lutions in the class is proved. 

1 General solution of evolutionary equation in 
qh; the form of power series and private cases. 

{3JQ| 1.1 Equations of spherical perturbations model 

As it is shown in Ref. [1] , spherically symmetric perturbations of Friedmann met- 
rics are described by one scalar function (5i/(r, rj) connected with perturbation 
of the metric tensor component 544 in the isotropic coordinates (r, rj) with the 
relation 



> 



6gn = a {r])6iy, (1) 
at that the representation turns out to be convenient 

^^^(r,ry)-/ifa)^ (2) 
ar ar 



- ^ , it allows to factor the solution into particlelike mode 

^ ' A,. ^^(^^) fo\ 
d- Sup = — , (3) 

equivalent to Newton potential of the point mass and the nonsingular mode 



^.0 = ^, (4) 
ar 

the "potential" of which ^(r, r/) and its first derivative by a radial variable at 
the beginning of the coordinates satisfy the relations 



lim |5'(r, 77)1 =]0, lim 



d'i'ir, 7]) 

r 



dr 



0. (5) 
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By a constant coefficient of barotrope k the mass of a particle source h{t]) and 
the nonsingular mode potential ^(r, 77) satisfy the evolutionary equations (75) 
and (84) [1]: 



2. ^^ + K)^^ 

77 (1 + 6K)^ rj^ 



77 (1 + 3k)2 772 ^ ' 

In the previous paper [1] an exact solution of the evolutionary equation ([5]) for 
the particlelike source mass was also obtained 

3(l + >i) 



^ = ^+771+3'= +M_77 1+3--, (1 + k)7^0; (8) 
fi^fi+V + Y' (l + '^) = 0. (9) 

In [1] it was noted that the peculiarity in the solution ([5]) by 1 + 3k = is coor- 
dinate disappearing by coming from the temporal coordinate i] to the physical 
time t 

„(1 + k) _ 1+3K 

Ja{ri)dr];^ t = cr;"^ i+3k ; ^ = (7^^3(1+^, (iq) 

Let us recall, the nonsingular part of the energy density perturbation Se(r, rj) 
and the radial velocity of medium v{r,r]) are determined with help of the po- 
tential function 4'(r, 77) by the relations 



Se 1 



So Anra^eo 



3-$-*"), (11) 



1 a $ 

^ ~^ Anra^eo dr r ^ ^ ^ 

where eo(^) is a nonperturbed energy density of Friedmann Universe? 

6(1+k) 2 6ft: 

£q ^ rj i+3k ; a '^771+31; Eoa^ ^ 7] i+3k . (13) 

In [1] by the variables separation method a general solution of the evolution- 
ary equation ([7]) satisfying the conditions ([5]) in form of an integral from Bessel 
functions was also obtained. In this paper we shall obtain more convenient 
solutions in form of power series. 
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1.2 Class C°° general solution in the perturbation area; 

K ^ 0, 1 + 

As it was shown in the previous paper [1], the singular part of the potential 
function Siy(r, rj) is uniquely extracted by the representation ^ in which the 
potential function 5'(r, t;) is nonsingular at the beginning of the coordinates, 
that is satisfies the relations ([S]) in consequence of which, in particular 

*(0,77) = 0. (14) 

Further supposing in the final neighborhood r — 0;r G [0,ro) in which the 
metrics perturbation is localized, the potential function 5'(r, 77) belongs to the 
C°° class, let us represent the solution of the evolutionary equation ^ satisfying 
the conditions ([5|) in the form of power series of the radial variable r 

00 

*(^,^) = I] (15) 



Let us underline the expansion (|15|) does not consist of a member with a zero 
power rin consequence of the relation (|14l) . Substituting the function \E'(r, ry)) 
in the form of the Eq. p3|) into the evolutionary equation ^ and setting the 
coefficients at each power rin the obtained equation equal to zero, we get a chain 
of linking equations 



* -n. * .."^^rn+l 6(1 + ^) ^2m+l 

- 0, ^2rn+l +^^^ (iTW 



= K(2m + 3)(2TO + 2)^'2™+3; m = 0,00, (k 7^ 0,-1). (16) 

Thus, by K 7^ the general solution of the evolutionary equation for the 
potential function ^(r, rj) is a series by odd powers of the radial variable r, that 
is by K 7^ 0, —1 the potential function \I'(r, ry) of the C°° class is an odd function 
of the radial variable r 

00 

M/(r,r;)=^M'2p+i(?7)r2?'+i. (17) 

p=0 

For the private particle solutions responding the specific physical conditions 
the series can be broken at any odd n ~ N ^ 3. In this case supposing for the 
last member of the series 

^'„(?7)=0; n>A^ = 2p+l, (p=l,2...), (18) 

we obtain from the Eq. (fT6)) a closed equation 



*,,„+2*S±l--«i±4^^.0^ (19) 

7/ (1 + 6K)'' T]'^ 
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As far as this equation does not differ from the evolutionary equation for the 
particleUke source mass at all, its general solution will coincide with the solution 
([HI) accurate within reterming 

2 3(1+k) 

^2p+l = C^7yl+3K + CP_^- 1+3K , (20) 

where and are some constants. 

Substituting this solution into the next to last equation of the chain p6p we 
shall get the equation for determining \I'2p-i: 



XT/ , o ^2p-l 6(1 + k) ^2p^l 
y2p-l +^ , Q \2 2 



K{2p + 3){2p + 2) 



2 3(1+k) 



(21) 



In consequence on this equation linearity its general solution is a sum of the 
general solution of the corresponding homogeneous equation 5'2p_i and a pri- 
vate solution of the inhomogeneous one ^'2p_i- But the general solution of the 
homogeneous one coincides with the already mentioned above solution (jS]) , and 
the private solution can be introduced in the form of a sum of two solutions 
corresponding to the two members of the right part p^ . Therefore, evidently 
the corresponding private solution has the form 

2 , ^ 3(1 + k) , ^ 



^-^p^i = Ap_i7/i+3-^^ + Bp_ir] 1+3- ^^ (22) 
Thus, we find 

Substituting the obtained solutions into the previous equations let us repeat 
the analogues calculations/ Thus we pointed out the algorithm of constructing 
a general solution of an evolutionary equation ((T]) reduced to repeating differ- 
entiation operations. This general solution corresponding to the highest power 
N = {2p+l) of the radial variable consists of 2N arbitrary constants appearing 
every time by solving the corresponding homogeneous differential equations. 

1.3 Case N=3 

Let us demonstrate the task solution in the simplest but as it turns out the most 
important case when N — 3 {p = 1) and the series ^T5\i consists of only two 
members corresponding to the valuesp — 0, 1. In this case from Eqs. (IT5]) - (l23l) 
we find 

2 3(1 + k) 

^-3 = Clryi+SK +CI77 i+3« ; (24) 
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2 3(1 + k) 



3.(1 + 3.)^,^.(|±M _ + ^^^^ 



1-3k 

2(7 + 9k) 2(1 -k) 



Let us study now the specific cases k — and 1 + k = which are out of the 
general solution (fT5|) . 

1.4 Nonrelativistic matter k = 

In this case we get from ([8]) the law of mass evolution of a particlelike source 
(see Ref. [4]) 

^ = + C_?7-^ (26) 

The equation ([7]) for the potential function ^ takes the form 

XE- + 2- - 6* = 0, K = 0, (27) 
whence we find 

^ = Mr)v^ + ^-ir)v'^ (28) 

where (/'±(r) are arbitrary functions r. In this case the mass of a particlelike 
source evolves according to the law 

M = M+7?' + M-'?"'- (29) 

1.5 Inflationary case k + 1 = 

In this case the equation for the field function F is elliptic 

l' + 2$- + $" = 0, (k + 1 = 0), (30) 
V 

and the radial velocity of perturbation is not determined by the equation (jl2p 
which in this case gives 

1(7)^°' i + -<'^ Pi) 

Integrating ([51]) we find: 

$ = (/.(r) + e(?7)r, (32) 

where 4>{r) and ^(77) are some arbitrary functions of their arguments. Substitut- 
ing this solution into the equation pip and dividing the variables we find the 
function $(r) 

<f = Ci - + r (^^2 + ^ + Ci77^ , (33) 

where Ci, C2 and C3 are arbitrary constants. 
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2 Retarding spherical perturbations in ultrarel- 
ativistic universe 

2.1 Boundary conditions for the retarding solutions 

By K > The spherically symmetric perturbations of Friedmann metrics are 
described by the hyperbolic equation, by k < - by elliptic one. Peculiarities 
of the hyperbolic equations ([7]) are convergent and nonconvergent waves 

r =F ^/kt] — Const, (34) 

spreading at the sound velocity 

Vs — ^/k- (35) 

Let us study the task for spherically symmetric solutions of linearized Einstein 
equations against the background of Friedmann metrics with zero boarding con- 
ditions for the potential function $(r, 77) at the sound horizon corresponding to 
the causality principle 

E: r = ro + \/k(77 - 7/o). (36) 

In terms of the introduced functions $(?', 77) and vE'(r, 77 these conditions can be 
written in the form 

In this case in consequence of the equations for perturbations out of the 
boundary of the sound horizon the perturbations of energy density, pressure 
and velocity must automatically turn into zerol^ 

An important private case of the boundary conditions ([57)) and ([55]) are the 
conditions at "the zero sound horizon " 

Eo : V^Tj. (40) 

In this case instead of (1571) and (155)) we have 

*('^.'?)I.v^, = mW; *V,^)I.v^, = o. (41) 



^ At the very boundary of the sound front the energy density perturbations and the veloc- 
ities may practically have final breaks. 
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2.2 Solutions with zero boundary conditions at zero sound 

Let us study some private solutions in the form of power series by radial variable 
satisfying the boundary conditions dH]). Such perturbations can be generated 
by the metrics fluctuations at zero moment of time and concentrated at this 
moment of time in zero volume. Coming to the most convenient radial variable 
Q 

r = y^P, (= -^P), (42) 

in the case of k = 1/3 let us write down the field equations ^ and ([7]). The 
dot denotes derivatieves by anew temporal variable 

/i + -/i - 4 = 0, 43 

2 • 2^* 

* + ^ - = 0. (44) 

then solving 

P = P+V + P-V'"^ , (45) 
where /z_ and are arbitrary constants. 

N=3 

As far as the Taylor-series expansion of the function ^'should not consist of the 
zero power r by definition and therefore all the even powers in this expansion 
automatically vanish, so setting further 

r/) = ^'i(77)r + *3(?7)r3 (46) 

and dividing the variables in the equation (|44p we get an equation for the 
functions5'i(?7): 

*i+f = 6*3; (47) 

^3+1^3-^ = 0- (48) 



In compliance with (j45p we find from 

^3^C+v + C^il-\ (49) 

where Cj. are some constants? Substituting the solution (|49| into the right part 
of the equation (j47|) and defining private solutions of the obtained equation we 
shall find its general solution 

vfi = C+rj + C^rr' + ^C+rj^ - 3C^ . (50) 
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Thus, 



(51) 



Further calculating ^'(p, 77)1^=,, = ^'(77, 77) , substituting the result into the 
boundary condition (|4ip and setting equal the coefficients in the obtained equa- 
tion by simultaneous powers 77 we shall get a set of equations for the constants 



_2 

77 





= A*-; 




cr 


= 0; 


77 




= M+; 


2 

77 


ct ' 


= 0; 


4 

7? 




= 0. 



(52) 



Thus, there are only two nonzero constants included into the expression for v. 
and — — l/2/i+. So, finally 



(53) 



4'(p,77) 



It is the obtained previously retarding solution with zero boundary conditions 
at zero sound horizon [2] andx(z) is the Heavyside function. At that supposing 
the following for the ultrarelativistic equation of state 



0(77)^77; (k=1/3), 
we obtain 



(54) 



3 -2 — 

rj rj-^ p 



xiv- p), 



(55) 



It is easy to see that in this case at the surface of zero sound front p — t] the 
boundary conditions (|4T1) are fulfilled identically. 



N=5 

Now let us study the fifth degree multinomial as a solution. In this case instead 
of the relations (H^ and (H5|) we have 

«'(p, 77) = ^-1(77)^ + *3(f/)p^ + 4'5(?7)P^ (56) 
^-i+f = 6^-3; (57) 

^-3+1*3 -21^ = 20^-5; (58) 
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(59) 



Similarly to the previous case we have 



(60) 



Substituting (|60p into the right part of the equation ([58]) , solving it similarly 
to the previous case and substituting the obtained solution into the equation 
(|57|) we get finally 



^3 = C+r, + C^r^-' 



2C+rj^ - IOC, 



C+rj^ - 3C3- 



(61) 
(62) 



and thus 

^(p, ry) = (C+77 + C^V^' + - + ^^5+7?'^ - I5C5- r;2)r+ 

+iC+v + C^ij-' + '^C+rf - IQC^y + [C+r^ + C^r,-^y. (63) 

Substituting the equation into the boundary conditions (PT|) and equat- 
ing the coefficients under equal 77 degrees we get equations for the constants 



(64) 



rj 





= M- 






= 0; 


V 




= 






= 0; 




-24C5- 


= 0; 






= 0; 




24^+ 
T^5 


= 0. 



Thus, 



0, 



(65) 



and the rest constants values coincide with the obtained ones above. It means 
that the solution coincides with the obtained one above. It is easy to show, 
adding of any new members of the series does not change the situation. 
Thus, we have proved the theorem: 



Theorem 1. The only spherically-symmetric solution of the class of 
linearized around Friedmann space-plane solution to the Einstein equations for 
an ideal ultrarelativistic fluid, corresponding to the zero boundary conditions 
at the zero sound horizon (HTj) . is tiie solution (|53p (it is equivalent to (55)). 
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2.3 Investigation of the retarding solution 



From the formula ([55)) one can immediately see continuity of not only the first 
radial derivatives but also the first temporal ones of the Sv metrics perturbation 
and the potential functions 



ddiy{p, ri) 



dp 



= 0; 



drj 



= 0; 



(66) 



(67) 



In Fig.l and 2 graphs of temporal evolution of the potential functions ^{p,i]) 
and —5v[p, 77) are shown 




Fig.l. Evolution of the potential 
function — $(p, 77). Bottom- up -q = 
1;2;3;4. Along the abscissa axis the 
radial variable values p = r^/i are 
laid off. 



Fig 2. Evolution of the metrics per- 
turbation —5v{p,rii). Bottom-up 77 = 
1; 2; 3; 4. Along the abscissa axis the 
radial variable values p = rVS are 
laid off. 



The second radial derivatives of the potential functions and metrics have a 
final break at the sound horizon. In consequence of this the energy density 
perturbation has a final break at the sound horizon also. Calculating the relative 
energy density of the spherical perturbations according to the formula ()lip we 
find 



£0 



1 



47r/9 \ri^ rj 
The jump of the relative energy density at the sound horizon is 



A 



£0 



(68) 



(69) 
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and it decreases by the time (see Fig. 3). Let us turn to the formula for 
the radial velocity (|12p . Substituting the equation for $ into this formula and 
coming to the radial variable p we obtain 



V 



1 + 2^ . (70) 



Evolution of the perturbation radial velocity is shown in Fig. 4. 



0.4 

Oi 






1 \ \ 


\\ 


P 

3 4 5 









"^-^ i 1 


-0.2 
-0.4 


1 

1 








\ 

\ 






-0.8 
-1.2 


\ 

\ 

\ 

\ 





Fig. 3. Evolution of the relative per- 
turbation of the energy density 5e/eo- 
From left to right 77 = 1;2;3;4;5. 
Along the abscissa axis the radial 
variable values p — r\/3 are laid off. 
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Fig 4. Evolution of the given ra- 
dial velocity v{p, ri)/ From left to 
right 7] = 1;2;3;4;5. Along the ab- 
scissa axis the radial variable values 
p = rVSare laid off. 



3 Conclusion 

Summarizing let us point out the following. It is easy to see that the obtained 
general retarding solution with a central singular source coincides with the ear- 
lier obtained ones [2], [3], [4] and [5]. However, note that in the mentioned works 
the solutions ([55)) - ((55|) were pointed out as private retarding solutions. The 
obtained retarding solutions belong to the class , however the second deriva- 
tives by radial variable have a break of the first genus at the sound horizon. 
As a result the energy density break at the sound horizon corresponds to the 
obtained retarding solution. Because of the pointed out conditions to find out 
the nature of the density break it is necessary to solve two tasks: to find the 
retarding solution for the barotrope arbitrary coefl[icient and solve the Cauchy 
task for the initially localized perturbation. We shall solve these tasks in our 
next papers. 
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